One-dimensional tight-binding models with correlated 
diagonal and off-diagonal disorder 
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We study localization properties of electronic states in one-dimensional lattices with nearest- 
neighbour interaction. Both the site energies and the hopping amplitudes are supposed to be of 
arbitrary form. A few cases are considered in details. We discuss first the case in which both 
the diagonal potential and the fluctuating part of the hopping amplitudes are small. In this case 
we derive a general analytical expression for the localization length, which depends on the pair 
correlators of the diagonal and off-diagonal matrix elements. The second case we investigate is that 
of strong uncorrelated disorder, for which approximate analytical estimates are given and compared 
with numerical data. Finally, we study the model with short-range correlations which constitutes 
an extension with off-diagonal disorder of the random dimer model. 
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I. INTRODUCTION 



The studies on the localization of electronic states in 
disordered solids can be traced back to the seminal paper 
of P. W. Anderson Q where the localization of quantum 
states was first discussed in connection with transport 
properties of a random lattice. This is a natural associ- 
ation, since the extended or localized nature of the one- 
particle wave-functions in a disordered system plays a 
key role in the determination of the metallic or insulating 
character of the system itself. As an outstanding exam- 
ple of the link between transport properties and wave- 
function localization in random media one could cite, 
for instance, the metal-insulator transition that occurs 
in disordered three-dimensional samples upon increasing 
the randomness of the material. 

The importance of the Anderson localization, however, 
far exceeds the limited field of condensed matter physics. 
That is mainly due to the fact that quantum interference 
is the fundamental mechanism responsible for the elec- 
tron localization in a random lattice; as a consequence, 
analogous localization effects can take place in any phe- 
nomenon that involves propagation of waves in a disor- 
dered medium Q . The same conceptual framework thus 
allows one to understand seemingly heterogeneous phe- 
nomena as the localization of water waves ||, the co- 
herent backscattering of photons j|, and the universal 
conductance fluctuations typical of mesoscopic samples 
(see, for instance, reference ||||). 

The basic one-dimensional (ID) Anderson model has 
proved to be a valuable tool to gain insight in the com- 
plex phenomenon of localization. The model is appealing 
because it is extremely simple and yet retains the capac- 
ity to provide a non-trivial description of the localization 
process. Mathematically, the Anderson model is defined 
by the tridiagonal Hamiltonian 



H 



\n)(n-l\ + \n)(n + l\ + \n)8 n (n\ 



(1) 



and disorder is introduced through the site energies 5 n 
which are supposed to be random and independent vari- 
ables (with an appropriate distribution). 

In spite of its undiscussed usefulness, the basic 
model (||) cannot account for many essential aspects of 
localization phenomena, such as the metal-insulator tran- 
sition. The twofold reason of this shortcoming lies in the 
ID nature of the system (Q) and in the statistical indepen- 
dence of the random site energies S n . As a consequence 
of these two key features, all eigenf unctions of the Hamil- 
tonian (Q) turn out to be localized, even if the disorder is 
infinitesimally small. 

The key role played by dimensionality in the localiza- 
tion processes led at first to the conclusion that only a 3D 
generalization of the Anderson model ([!]) could reproduce 
such basic characteristics of disordered solids as the exis- 
tence of mobility edges or the metal-insulator transition. 
The progress in the investigation of ID systems, however, 
eventually led to the realization that such models can ex- 
hibit a richer behaviour than it was previously thought. 
In fact, it was discovered that even in one-dimensional 
lattices there are classes of random potentials which al- 
low for extended states (see, e. g., reference 0). In all 
these ID systems the disorder exhibits spatial correla- 
tions, in contrast to the totally uncorrelated potential of 
the standard Anderson model. A significant example of 
a system where the correlations of the site potential can 
produce delocalized states is given by the so-called 'ran- 
dom dimer', which is characterized by a peculiar form of 
short-range correlations of the potential ||. 

The analysis of models with correlated disorder revived 
the interest in ID chains and revealed the essential role 
of correlations. One must observe, however, that the 
short-range correlations of the random dimer and sim- 
ilar systems can give rise at best to a discrete set of 
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extended quantum states: in other words, the localiza- 
tion length diverges only for discrete values of the energy. 
The situation may change when long-range correlations 
come into play, as it was recently discovered in J^] , where 
the general case of diagonal disorder with arbitrary cor- 
relations was considered and a direct relation between 
the localization length and the potential pair correlators 
was established. Using this relation it was shown how 
to reconstruct a site-potential giving rise to any given 
form of energy dependence of the localization length || . 
A particular conclusion is that even ID random lattices 
can possess a continuum of extended states (and mobility 
edges), provided the disorder exhibits appropriate long- 
range correlations. 

In the present work we follow an approach similar to 
that of reference Q by extending the study to a general 
case of ID Anderson model with any kind of diagonal 
disorder and off-diagonal nearest-neighbour interaction. 

This work is organized as follows. The rest of the intro- 
duction is devoted to the definition of the model under 
study. In Section || we analyse the case in which the 
site-energies and the site-dependent part of the hopping 
amplitudes are small. Although we focus our discussion 
on the case of random site- and hopping energies, the 
validity of our results extends to any non-random case 
as well. In Section III we consider the case of strong 
and uncorrelated disorder, deriving an approximate ex- 
pression for the localization length and comparing it with 
the numerical data. In Section IV we investigate the case 
(reminiscent of the random dimcr model) in which the 
disorder (both diagonal and off-diagonal) has arbitrary 
intensity and exhibits peculiar short-range correlations. 
The concluding remarks are exposed in Section [v| 



A. Definition of the model 

We study the localization properties of the following 
one-dimensional tight-binding Hamiltonian 

H = Y^ [l«>7n-i<n - 1| + \n)j n+1 (n+l\ 

n 

+ \n)S n (n\] (2) 

where the site energies S n and the off-diagonal elements 
7„ are arbitrary real variables. From the physical point 
of view, the Hamiltonian (||) describes an 'electron' mov- 
ing in a discrete ID lattice. The energies S n measure the 
strength of the bond of the electron to the 'atoms' of the 
chain, while the hopping terms 7„ represent the proba- 
bility amplitudes for an electron localized in a single site 
to jump to the nearest atoms to the right or to the left. 
Notice that in our model the transition amplitudes 

A n ^ n ±i = (n ± l\H\n) = j n 

depend only on the initial state of the electron and not 
on the final one. This physical feature is connected to 



the non-Hermitian character of the operator (g). With 
regard to this point, we stress that, although the oper- 
ator (^]) is not Hermitian, its eigenvalues are real: this 
allows one to consider it as a physically sound Hamilto- 
nian. 

In order to show that the eigenvalues of H are real, 
let E be an eigenvalue of the operator (||) and \iPe) the 
corresponding eigenvector 



H\rl> E ) =E\ip E ] 



(3) 



Let us now consider the operator H T , i.e. the transpose 
of H and let us denote with \4>e) the eigenvector of H T 
to the eigenvalue E, so that one has 



H T \4>e ) = E \4>e)- 



(4) 



Notice that if E belongs to the spectrum of H, then the 
same value is also an eigenvalue of H T , since a matrix 
and its transpose share the same spectrum. Multiplying 
equation (Q) by (4>e\ one can obtain 



E(1> E \<t>E) = (iPe\H t \cI>e) = 

(i^l <M = #*(^M <M 



(5) 



where we have used equation (^) and the fact that if is a 
real matrix (and therefore its transpose coincides with its 
adjoint). Equation (|J) implies that E = E* , unless one 
has (ipE I 4>e) = 0. One can rule this last possibility out, 
however, due to the specific structure of the operator (||). 
Indeed, for the Hamiltonian (0) one has that 



H T = A HA 



(0) 



Relation (g) between H and 
1 \iPe), and this allows one to 



where A = J2 n \ n )^( n \- 
H T implies that \4> E ) = A 
exclude the possibility that \i/je) and \cj) E ) are orthogo- 
nal. This leads to the conclusion that E = E* , i.e. that 
the Hamiltonian (|^) has real eigenvalues. 

Below, we consider the stationary Schrodinger equa- 
tion H \ip) = E with eigenvectors 



|^> = ^ n |n). 



Using the explicit form (g) of the Hamiltonian, one can 
easily see that the Schrodinger equation for the ampli- 
tudes ip n takes the form 



7n+l^n+l + 7„_l-0n-l = (E - S n ) tj) n . 
It is convenient to introduce the new variables 



and 



4>n = lni>n 

E-S n 



£n (E) = 



In 



(7) 



(8) 



(9) 
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so that equation (|7]) can be cast in the simpler form 

4>n+l + <f>n-l = £,n4>n- (10) 

The introduction of the amplitudes (^) and of the site- 
energies (^) thus formally allows one to reduce the prob- 
lem of determining the eigenstates of model (|J) to the 
equivalent problem of studying the zero-energy eigen- 
state of the tridiagonal model ( |l0| ) with diagonal-only 
disorder. According to this interpretation, for any en- 
ergy value E in the original Schrodinger equation there 
is a corresponding realization {£„ (E)} in (|l0|). 

Having thus defined our problem, we proceed to solve 
it for three distinct physical cases: the limiting case of 
weak disorder (with arbitrary correlations), the opposite 
case of strong (uncorrelated) disorder and the case of dis- 
order of arbitrary strength with short-range correlations 
('iV-mer model'). In the following, we will consider off- 
diagonal terms of the form 



In 



1 



(11) 



The reason for representing the energies j n in the 
form (|ll] ) is that in this way we separate two physi- 
cally distinct contributions to the off-diagonal matrix el- 
ements: the first one originates from the kinetic energy 
term of the Hamiltonian, while the second one, e n , repre- 
sents the fluctuating part of the interaction between near- 
est neighbours. The latter contribution is the source of 
randomness for the hopping energies j n of our model (||). 



II. WEAK CORRELATED DISORDER 



A. Localization length 



In this Section we consider the general case of weak 
disorder: 



|e„| < 1 and \S n \ < 1 



(12) 



for any kind of site-energies S n and interaction e„. 

Notice that, since the randomness of the hopping terms 
7„ comes only from the interaction energies 6 n , the con- 
dition of weak disorder for the off-diagonal terms is given 
by the first of the relations ([[2]) (and not by the condition 
7„ < 1). Under the conditions ( |l2| ) one has 

£„ ~ E - (Ee n + S n ) 

so that equation (|l^) can be approximated by 

4>n+l + 4>n-l = (E - Ee n - S n ) <f) n . (13) 



1. Classical representation of the quantum model 

An effective way to study the quantum model dl3| ) 
consists in representing it in terms of the classical two- 
dimensional Hamiltonian map 



Pn+l 
Xn+l 



with 



and 



(p n + A n X n ) COSfi 4 

- (p n + A n x n )smfjL- 



E = 2 cos \i 



x n sin fi 
- x n cos jJL 



A n 



S n + 2e„ cos \i 
sin fi 



(14) 



(15) 



(16) 



The map (14) describes the behaviour of a harmonic os- 
cillator subjected to periodic delta kicks of amplitude 
A n . The map can in fact be derived by integrating over 
a period T — 1 the equations of motion of the kicked 
oscillator defined by the Hamiltonian 



H 



A n S(t-n). 



The complete equivalence of the quantum system (|l^ ) 
with the classical one defined by the map (|IJ) can be 
easily proved by eliminating the p variable from equa- 
tions Jl4|). In this way, one obtains the relation 

X n+ 1 + X n -l = (2 COS fi - A n SUl fl) X n 

which coincides with equation (|l^) if one identifies the 
variable x n with the site amplitude 4> n and if the param- 
eters of the quantum model (13|) and those of the classi- 
cal map (|l4|) are related by the equations Jl5| ) and Jl6|). 
Therefore it becomes possible to analyse the solutions of 
equation ( |l3| ) in terms of 'trajectories' in the phase-space 
of the map (|TJ). In such an approach, localized quan- 
tum states correspond to unbounded trajectories in the 
classical phase-space, while extended quantum states are 
represented by bounded trajectories (see details in JTo||). 

It is convenient to express the map (|lj) in polar coordi- 
nates (r, 0) introduced via standard relations p = r cos 
and x = rsmO. The substitution in equation ( |l4| ) gives 

9 n - A*) + An sin 9 n cos fj] 
9 n - fi) - An sin 6 n sin /x] 



sm6> n+ i = jy- [sin 



l 

D„ 
1 

D„ 



where 



Dn 



r n +i 



1 + A n sin (20 n ) + A% sin 2 



Using this approach, the inverse localization length I 1 
(or Lyapunov exponent A) can be expressed as 



1 N 

r 1 = A = lim — Vlo 

AT^oo N ^ 



Xn+1 



lim — los 



JV 



N^oo N 



r n +i sintf„ +1 
r„ sin n 



(17) 
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Except that at the band edge (i.e. for \E\ — > 2), expres- 
sion (flTh can be safely reduced to 



ao = E* e " 2V( ^ 1) - 



(23) 



k=i 



A = lim - Y log ( ) = lim 1 V \ogD n (18) Substituting @ into @, one can write 



N->oc N 



(see Q for details). 

For weak disorder, the logarithm in (|l^) can be ex- 
panded and to the second order of perturbation theory 
one gets 



A=±<^) + ^ n siii(20 n )> 



(19) 



where the symbol (• ■ •) stands for the 'time'-average, 

~N 



2. Computation of the noise-angle correlator 

To compute the correlator (A n sin (26 n )) with second- 
order accuracy we follow the approach of || . We define 
the noise-angle a k and the noise-noise q k correlators as 
follows, 



2i 



ak 



2i V/A „2i9„_ fc \ 



-^-e 2lfl {A 



(20) 



and 



Qk — ~ o~\A n A n — f-) 

_ E 2 (e n e n _ k ) + 2E(e n S n _ k ) + (6 n 6 n _ k ) 
where o~ 2 A is the 'noise'-variance 



(21) 



2v E 2 (e 2 n )+2E(e n S n ) + (S 2 n ) 



1-E 2 /A 



The terms (g0j) measure the 'temporal' correlations of 
the noise A n with the angle variable ni while equa- 
tion ( ^l| ) defines the normalized autocorrelation function 
for the noise variable (|l6|). Notice that, since the noise 
strength (|l^) is a function both of random terms e n and 
S n and of the energy E = 2 cos fi, the noise-noise cor- 
relator (^lj) depends on the energy E as well as on the 
(spatial) correlators of the random variables S n and e n . 
From the definitions (|2^), (|2l]) it follows that 



(A n sin(20„))=Re( -fe~ 2 ^a 



(22) 



As in reference M it can be shown that the correlators 
a k satisfy the relations 

Ofc-i = e~ 2lfJ -a k + q k , for k = 1, 2, 3, . . . 

which in turn imply that 



(A n sin (20„)) = ^ ]T q k cos (2/xfc) 
fc=i 

so that the Lyapunov exponent ( |l9| ) takes the form 

A = A ^(/x) (24) 

where 

Xo{E)= < = m±^M^m (25) 



and 



8(1 - £ 2 /4) 



y(/i) = l + 2^^ (/i) cos (2//fc) 

k=l 



(26) 



while the parameter fi is related to the energy E through 
equation (iq). 



3. Discussion of the results 



Expression (|24j) , together with equations (|25| ) and (|2f 
gives the localization length as a function of the energy 
and of the diagonal and off-diagonal potential correla- 
tors. Notice that this is a very general result, because its 
validity rests exclusively on the weak disorder assump- 
tion ( p^ ) and does not depend on the particular form of 
the variables 6 n and e n . Equation (|^) can therefore be 
applied to a broad variety of specific problems. 

A simple examination of formula (|24|) also reveals that 
the inverse localization length A is the product of two fac- 
tors, An and (p (fi), defined respectively by equation ( g^ ) 
and (|26j). This factorization is physically meaningful, be- 
cause An represents the Lyapunov exponent for the case 
of totally uncorrelated disorder (i.e. when q k = for 
k > 1), while the function tp (/x) describes the correction 
introduced by the spatial correlations both in diagonal 
and off-diagonal terms. 

Let us now analyse separately the two factors fl25| ) 
and (^6|) . As was indicated, the first one gives the inverse 
localization length when no correlations exist among the 
matrix elements of the Hamiltonian (||). This system 
corresponds to a simple generalization of the Ander- 
son model, in which not only the site energies but also 
nearest-neighbour hopping energies are independent ran- 
dom variables. From this point of view, the expres- 
sion ( p5[ ) can be seen as an extension of the known for- 
mula for the standard Anderson model, 



A = 



W 



96(1- £ 2 /4)' 
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Indeed the inverse localization length (|25|) reduces to this 
form if one assumes that there is no off-diagonal disorder 
(i.e. (e^) = and (e n 5 n ) = 0) and that the variance of 
the site energies is the one fixed by the Anderson distri- 
bution (i.e. (S n ) = W 2 /12). Expression (p5[), however, 
also shows that the introduction of off-diagonal random- 
ness modifies in a non-trivial way the energy dependence 
of the Lyapunov exponent with respect to the standard 
diagonal case. In particular, extended states can arise for 
specific values of the energy: if, for instance, the diago- 
nal and off-diagonal disorder are identical (e„ = 8 n ), the 
localization length diverges for E = — 1. We discuss this 
feature in more detail in the next Subsection. 

The factor ( p6| ) accounts for the modifications gener- 
ated by the correlations both in the diagonal and off- 
diagonal terms. One should stress that this factor de- 
pends both on the correlations and energy. This is a key 
feature, because it implies that, as an effect of the corre- 
lations, the sample can become more or less transparent 
for specific values of the energy. As can be seen from 
equation if tp (fi) > 1, the Lyapunov exponent is 

increased with respect to the uncorrelated case, mean- 
ing that correlations make the sample more opaque. On 
the contrary, when tp(fi) < 1, the chain becomes more 
transparent. Thus, if the function (|6|) goes to zero in a 
certain energy range, the states of the system will be ex- 
tended inside that region and localized outside of it: this 
shows how mobility edges can arise even in ID systems, 
provided that the potential has the proper correlations. 

Once the statistical properties of the potential are 
known, formulae (|2^), fl25| ) and (26) exactly determine 
the Lyapunov exponent as a function of the energy. They 
can be used, however, also to solve the inverse problem, 
i.e. to specify what the correlations of the disorder must 
be for the Lyapunov exponent to exhibit a given energy 
dependence. This is a problem of practical interest be- 
cause of the relevance of its potential technological ap- 
plications: knowing what particular disorder generates 
a specific energy dependence of the Lyapunov exponent 
might lead to the construction of superlattices with the 
required transport properties. 

Let us now write explicit relations between statistical 
properties of the disorder and the Lyapunov exponent for 
the case when the diagonal disorder is independent from 
the off-diagonal one ( (e n S n -k) = ). For this we write 
the function (p (/z) in Fourier series, 



+oo 



(m) = Lpk e 



i2fik 



(27) 



k= — oo 



where the expansion coefficients are defined by the rela- 
tion 



1 f~ 



9k 



if (m) e- t2f * k d/i 



Inserting ( p7| ) and ( pl| ) in equation (26) one obtains the 
relations 



{5 n )fk + (e n ) (<Pk-i + 2ip k + (fk+i) = {S n S n -k) 

) (28) 

which constitute the explicit link between the Fourier 
components of ip (/i) and the disorder correlators 
{6 n S n ~k} and (e n e n ~k)- Once the Lyapunov exponent 
(and hence the factor tp (//)) is known, the relations ( pgj ) 
set a constraint on the statistical properties of S n and 
e„. Notice that, in contrast to the purely diagonal case 
discussed in [||, here the interplay of diagonal and off- 
diagonal disorder does not allow a complete specification 
of the correlators from the knowledge of the localization 
length. 

The relations ( p8|) admit a particularly simple solution 
for the case in which the diagonal and off-diagonal disor- 
der are statistically independent (i.e. (e„S n ^k) = 0), but 
share the same statistical properties (so that their au- 
tocorrelation functions are equal: (5 n 5 n -k) = (tntn-k))- 
In this peculiar case, the relations ( |28|) do determine the 
disorder correlators: indeed, one has 



(SnS n -k) 

(SI) 



(4> 



fk- 



This result can be understood from the observation that 
in this special case the noise-noise correlators ( fH| ) take 
the form = (S n d n -k) / and no longer depend on 
the energy, so that they can be interpreted as Fourier 
coefficients of the series (26). 



B. Numerical data 

We checked numerically the validity of (J|J) using the 
classical map approach for the computation of the Lya- 
punov exponent. We consider the map 



Pn+l 

X n +i 



Mr. 



with 



Mn = 



COS /! 



A n cos [i + 7„ sin ji 



~ : i^ sinfJ ' 7^ (in cos fi ~ A n sin ^) 



(29) 



(30) 



which is exactly equivalent to the original Schrodinger 
equation (H), provided the parameters E and A n of the 
matrix (|30| ) are defined by the relations (|l5| ) and (|l(|). 
This means that, if ( |l5| ) and ( |l6| ) are fulfilled, then the 
values of the X coordinate at the time steps n — 1, n 
and n + 1 obey the same relation expressed by (^) for the 
probability amplitudes ip n . 

It should be pointed out that the determinant of the 
single time-step map ( |30| ) is not unitary, detM„ = 
7n/7n+i and, therefore, the determinant oscillates 
around the unit value as a function of n. This does not 



mean, however, that the map (29) does not conserve the 
total flux. Indeed, the disordered sample represented by 
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the Hamiltonian (g) cannot actually be infinite; therefore 
some kind of boundary conditions have to be imposed. 
A typical boundary condition is the periodic one: in this 
case one assumes that the disordered sample is a closed 
chain made up of N 'atoms' so that the site- and hopping- 
energies obey the equations: S^ + i = Si and 7at+i = 71. 
Consequently, one has 



TTdetM^Ii.lH 

72 73 



k=l 



In-1 In _ 71 



(31) 



so that the total determinant across the sample is ex- 
actly equal to one. A more physical boundary condition 
arises if the disordered sample is embedded between two 
perfect leads. This physical condition translates in the 
mathematical requirements 8 n — and 7„ = 1 for n > N 
or n < 1, which, in turn, lead again to the condition j3l| ) 
for the total determinant of the map through the whole 
sample. 

Let us now see how the map (|2^) can be effectively 
used for the computation of the Lyapunov exponent. As 
was shown above, it is convenient to introduce the po- 
lar coordinates X — r sin 9 and Y = r cos 9 and cast the 
map (p9|) in the form 



cos 9 n+ i = [cos (9 n - fx) 

+ (A n cos /1 + e n sin //) sin 9 n ] 
sin 9 n+1 = ^ I sin ( e n - M) 

+ (A n sin \i + e n cos fi) sin 9 n ] 



(32) 



which allows an easy computation of the expression Jl8| ) 
for the Lyapunov exponent. In this way, we tested the 
validity of the formula (^J) for the special case in which 
the hopping energies e„ are independent random variables 
with a common distribution, while the site-energies are 
defined by the relation 



S n = ae n with lal < 2. 



(33) 



Taking into account Eqs. (|33j) and ([11]), the Schrodinger 
equation (0) can be written as 



(1 + e n+ i) ip n+ i + (1 + e„_i) tpn-i = (E - ae n ) tp n . 



The box distribution 



(34) 



(35) 



of width W was chosen for the variables e n . For this par- 
ticular model the general expression (^) reduces to the 
following form; 



A = 



W 2 (E + af 
96 (l-E 2 /4)' 



(36) 



This formula, among other things, shows that for the 
above-defined model the state corresponding to the en- 
ergy value E = — a is always extended, regardless of the 
noise strength fixed by the parameter W . In particu- 
lar, when a = (i.e. there is no diagonal disorder), the 
model ([m]) exhibits an extended state at the band centre, 
i.e. for E = 0. The value a = 1 corresponds to the case, 
cited in the previous Subsection, in which the binding 
and hopping energies are equal (e„ = S n ) and one has 
a transparent state for E = — 1. The comparison with 
numerical data revealed an excellent agreement between 
the theoretical prediction ( |36| ) and the actual behaviour 
of the Lyapunov exponent A, as can be seen from Fig. |l| 
and Fig. || which represent A as a function of the energy 
E for two values of the parameter a. 



Numerical data 
Theoret. expression 



0.008 - 



0.006 - 



0.004 - 



0.002 - 




0.5 1 1.5 2 

Energy 

FIG. 1. Lyapunov exponent vs. energy. The figure refers 
to the case a = and W = 0.1. 
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FIG. 2. Lyapunov exponent vs. energy. The figure refers 
to the case a = 1 and W = 0.1 
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III. STRONG (UNCORRELATED) DISORDER 
A. Analytical estimates for the localization length 

For the strong disorder case we consider the model (|^) 
with diagonal terms of the form (|33|), i.e. S n — ae n with 
\a\ < 2 . The Schrodinger equation is then given by ( p^ ) 
and the model is completely dehned once statistical prop- 
erties of the random variables e n are specified. Through- 
out this Section we assume that the energies e„ are in- 
dependent random variables with a common probability 
distribution given by the box distribution ( pq ) of width 
W. The case of strong disorder is then defined by the 
condition W ^> 1. 

Again, by introducing the rescaled amplitudes Mj and 
the binding energies ([}]), we can reduce the Schrodinger 
equation ( p4| ) to the form (|l^) with only diagonal disor- 
der. In the latter model the site-energies (O) take now 
the form 



By introducing polar coordinates through the relations 
4>n = r n cos9 n and 4> n -i — r n siii9 n , one can write the 
map (MOT) in the form 



cos 9 n+ i = sin 9 n 

sin 9 n+ i = [cos 9 n + sin 9 n ] 



(42) 



where 



D n = 



r n +i 



r.„ 



As before, the Lyapunov exponent is given by the rela- 
tion ( O ) and can, therefore, be expressed by the integral 



A = y dep(0 J d9 Pi (9) log 



(43) 



E 



1 



and are independent random variables whose distribution 
p (£„) is related to the probability p (e n ) of the original 
variables e„ by 



where the symbol (9) stands for the invariant measure 
of the variable 9 of the map ([42]), and the subscript £ 
stands in order to stress that such a measure depends on 
the noise distribution p(£ n ). 



E\ (E-i n 

—p 



in 



B. Decomposition of the invariant measure 



In explicit form, when p (e„) is the Anderson distribu- 
tion fl35|) , the corresponding probability for the variables 
£„ reads 



where 



and 



! - LL| 1 if ^ < or ^ < C„ 

if < Zn < ^ 



;■(!) 



2g + aiy 2£-aiy 
2-W ' 2 + W 



(37) 



(38) 



The integral 
1Li_ + J_ with 



31) can be separated in two parts, A 



I± = J diP± (0 / d»p± (0)log£>(£, 



Here p± (£) are the noise distributions 



P+ (0 = 



_/p(0 for |C| >2 







elsewhere 



i 



(2) 



max 



2£ + aTK 2£-aW 
2-W ' 2 + H/ 



(39) 



Following the Hamiltonian map approach already used 
in the previous Section, we can cast equation ( |l0| ) in the 
form of the two-dimensional map 



■>n+l 



where T„ is the transfer matrix 



^ -1 
1 



(40) 



(41) 



and 



v- (0 



P(0 for \i\<2 
elsewhere 



while p + (9) and p_ (0) are the corresponding invariant 
measures. This decomposition is meaningful because nu- 
merical computations have shown that, upon increasing 
the noise strength W, the invariant measure p_ (S) be- 
comes more and more sharply peaked around the single 
values 9 — ixh with k integer (as can be clearly seen from 
the Fig.i). 
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FIG. 3. Invariant measure p- (6) for various noise 
strengths (data refer to the energy value E=0.01 and a = 0) 



Once more, we make use of the fact that the 
Schrodinger equation is exactly equivalent to the 
2D map ( f40| ) whose eigenvalues are 

A (l,2) = (45) 

For |£ n | < 2 the eigenvalues ( p5| ) take the form A« = 
e ±*M« with £„ = 2 cos so that, at each step, the evolu- 
tion dictated by the map (^) results in a simple rotation 
in the two-dimensional phase-space. For > 2, how- 
ever, both eigenvalues of the map (^) are real and this 
consideration allows one to compute the Lyapunov expo- 
nent ( |44| ) using the approximated expression 

A=(log|A|)= /<^ + (0logA(0 (46) 



As a consequence, the integral /_ becomes negligible 
for strong disorder. The formula ( f43| ) can therefore be 
approximated by 



^ — / dCp(0 

iei>2 



d9 P+ (6)]ogD(£,6). 



(44) 



The analysis of this expression reveals the simple de- 
pendence A oc l/W of the Lyapunov exponent on the 
noise strength. In fact, in ( |44| ) the VF-dependence of 
the integral is introduced in a twofold way, through the 
functional l/W dependence of the noise distribution j37| ) 
and through the integration range for the £- variable. For 
large enough W, the integration range is simply deter- 
mined by the condition |£| > 2 because the limits (J3q) 
and ([39]) of the interval where the distribution (|37|) is 
zero are located inside the interval |£| < 2. Therefore, 
for the integral (44) the only dependence left on the noise 
strength is the one due to the factor 1 /W contained in the 
integrand. Physically, this means that, for strong disor- 
der, the localization length grows linearly with the noise 
intensity. In contrast to the standard Anderson model, 
therefore, an increase of the lattice disorder reduces the 
localization of the electronic wave-functions instead of 
enhancing it. 



C. Another approach for the computation of A 

Although formula (fht ) allows one to predict the rela- 
tion of inverse proportionality between A and W without 
carrying out any explicit calculations, the determination 
of the proportionality constant requires a direct evalua- 
tion of the integral. Unfortunately, the analytical com- 
putation of integral (Q) is out of reach, since there is 
no way to determine the invariant measure p+{0). One 
can overcome this difficulty using a different approach 
successfully used in |11 for an analogous problem. 



where A = max { I A (1) I , I A (2) I } . The evaluation of inte- 
gral (E3) leads to the formula 



A = 7 (a) 



\E + a\ 
W 



(47) 



where 



The validity of the formula ( [47[ ) was numerically 
checked; the computations showed a good agreement be- 
tween the theoretical prediction and the actual behaviour 
of the Lyapunov exponent, as can be seen, for instance, 
from Fig. ^ and Fig. || 
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FIG. 4. Lyapunov exponent vs. energy for various noise 
strengths with a = 
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FIG. 5. Lyapunov exponent vs. energy for various noise 
strengths with a — 1 



Notice that the Lyapunov exponent ( |47|) exhibits the 
expected 1/W dependence on the noise strength. As we 
explained, this is a consequence of the fact that the re- 
gion < £ < £( 2 ) where the integrand of formula fl46| ) 
is zero lies outside the integration domain |£| > 2. In 
passing we remark that, if the condition \a\ < 2 were 
not fulfilled, this conclusion would not be true. There- 
fore, for |a| > 2 the integral ( f46| ) would depend on W 
not only through the integrand but also through the lim- 
its of the integration domain. The computation of the 
integral remains possible also in this case, but then the 
result is no longer as simple as the one expressed by for- 
mula (47). That is why we chose to restrict the range of 
the parameter a to the interval \a\ < 2. 



IV. DISORDER WITH SHORT-RANGE 
CORRELATIONS 

In this Section we focus our attention to the case in 
which the diagonal and off-diagonal disorder in the gen- 
eral model (||) exhibit short-range correlations. More 
precisely, we consider the model defined by assuming 
that the couples of random variables ( n = (6 n , e n ) in the 
Schrodinger equation ([7|) fulfill the following conditions: 

1. The energies £„ can take only two values: £+ = 
(<5+,e+) or C- = (S-,e-). 

2. In every succession {£ n } the values and ap- 
pear in iV-uples made of N consecutive identical 
terms. 

3. The iV-uples (£+, ...,£+) and (£_, ...,£_) appear 
at random with equal frequency along the succes- 
sion {£„}. 

If N = 2, the pairs (£+,£+) and (£_,£_) are called 
dimers; for N — 3 the corresponding triplets are named 
trimers, while for generic N one speaks of 'iV-mers'. 



Making use of the amplitudes (||) and the energies (j|) 
as in the previous cases, one can write the Schrodinger 
equation (0) in the form where the diagonal ener- 
gies (||) now take two values 



E-S± 
l + e± 



(48) 



and appear in equiprobable iV-uples randomly positioned 
along the succession {£ n }. In other words, the model (||), 
with both diagonal and off-diagonal disorder specified by 
the preceding 1), 2), and 3) conditions, can be associated 
with the iV-mer model already studied in the literature 
(see Q and references quoted therein). 

As in reference |l(| , we can use the equivalence of the 
model ( |Io| ) with the map ( |4(i|) and study the solutions of 
the Schrodinger equation (|10[) in terms of the trajecto- 
ries of the map (j4Cf). With this interpretation in mind, 
we can address the problem of extended states for the 
model ( |lo| ) by considering a single iV-mer of type em- 
bedded in an infinite chain of N-mers. In this case the 
orbits of the map (|4(]) are not affected by the single 
iV-mer provided that the total transfer matrix through 
the iV-mer is equal to the identity matrix (apart from a 
plus/minus factor). In other words, one must have 



= ±E 



(49) 



where T + is the transfer matrix ( pll| ) (with £ n = £+) and 
E is the unit matrix. 

The 'transparency' condition ( [49| ) can be satisfied pro- 
vided that the stability condition |£ + | < 2 for the 
map (40) is fulfilled. In this case the eigenvalues of the 
matrix T + have the form e ±i ' 1 + (with £ + = 2cos/i + ) and 
the equation ( [49| ) is equivalent to exp(i(i + N) = ±1 which 
in turn leads to /i + = nk/N with k integer. On the other 
hand, the stability condition |£ + | < 2 requires to discard 
the k values k = m ■ N with m integer. Using (48) and 
the relation £ + = 2 cos /i+ , we come to the conclusion 
that the ./V-mer is transparent for the following energy 
values 



E 



2(1 



with k = 1,2, ... ,N - 1. 

The study of the single iV-mer case allows one to under- 
stand that, in the general case of a chain where A^-mers 
of type and randomly alternate, the model (||) with 
disorder correlated by A^-uples has extended states when 
the energy takes the values 



E 



5+ 
<5_ 



2(1 
2(1 



£_)cos(at 



(7vk\ 

77 k 



(50) 



where fc = 1, 2, . . . , JV — 1. In the particular case when 
the diagonal and off-diagonal disorder are related by the 
linear equation S n = ae ni one should add the special 
value E = —a to the list (BQ), because in this case the 
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variables (|48|) take the single value 
disappears from equation (|Io|). 



-a and disorder 



The theoretical prediction (50) was checked by nu- 
merically computing the Lyapunov exponent for the sys- 
tem (^) with dimer- and trimer-correlated disorder as 
a function of the energy. Numerical experiments have 
confirmed the validity of formula (|5(i|), specifically, the 
Lyapunov exponent vanishes whenever the energy takes 
one of the values j50|), see Figs. [| [7|. 
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FIG. 6. Lyapunov exponent for dimer-correlated disorder 
with 5± = 0, e+ = 0, e_ = 0.5 
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FIG. 7. Lyapunov exponent for trimer-correlated disorder 
with <5± = 0, e+ = 0, e- = 0.5 



main object of this study is the analysis of the localization 
length and its dependence on underlying correlations in 
the diagonal and off-diagonal potential. By making use 
of the reduction of the quantum model to a proper clas- 
sical Hamiltonian map, we have considered a few most 
interesting cases of the interplay between diagonal and 
off-diagonal correlated disorder. 

The first case which has been analysed in details is the 
model with a weak site-potential and a weak fluctuat- 
ing part in the off-diagonal matrix elements. In this case 
we have derived a general analytical expression for the 
localization length valid for any kind of randomness in 
the potential. More specifically, the localization length 
is expressed in terms of pair correlators for the diagonal 
and off-diagonal matrix elements. Using this expression, 
one can consider any kind of correlated disorder, as well 
as deterministic potential and hopping off-diagonal am- 
plitudes. The expression obtained for the localization 
length allows to reveal how this latter quantity is influ- 
enced by pair correlations. A few specific cases of corre- 
lated disorder have been discussed to illustrate the role 
of the correlations. One of the most interesting results 
is that if the diagonal and off-diagonal disorder are pro- 
portional to each other, there is a specific value of the 
energy for which the eigenstates are extended. For finite 
sample, therefore, there is a region where all states are 
fully transparent as in the dimer models with short-range 
correlations. 

Another case studied in this paper is that of strong 
disorder (both diagonal and off-diagonal). In this case 
we have been able to obtain an approximate analytical 
expression for the localization length which shows that, 
when lattice randomness increases, localization effects 
weaken. This feature contrasts with the behaviour of the 
standard ID Anderson model (with only diagonal disor- 
der) and is a non-trivial consequence of the interaction of 
the site-potential with the off-diagonal hopping energies. 

Very recently [ |12[ an analytical expression for the in- 
variant density p^(p) in the strong disorder case has been 
obtained in a different context. 

Finally we have examined the case of disorder with 
short-range correlation, focusing our attention on the 
specific N-mer model, where disorder correlations are in- 
troduced by assigning the same value of the site- and 
hopping-energies to N consecutive chain sites. For this 
case we have showed how correlations give rise to a dis- 
crete set of extended states. 
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